A^-particle scattering matrix for electrons interacting on a quantum dot 



OO 

o 
o 

(N 



CD 



o 

> 
OO 

o 



o 



X 



A.V. Lebedev°^^ G.B. Lesovik^ and G. Blatter ° 
Theoretische Physik, Schafmattstrasse 32, ETH-Zurich, CH-8093 Zurich, Switzerland and 
^L.D. Landau Institute for Theoretical Physics, RAS, 119334 Moscow, Russia 

(Dated: June 10, 2008) 

We present a non-perturbative expression for the scattering matrix of N particles interacting inside 
a quantum dot. Characterizing the dot by its resonances, we find a compact form for the scattering 
matrix in a real-time representation. We study the transmission probabilities and interaction- 
induced orbital entanglement of two electrons incident on the dot in a spin-singlet state. 

PACS numbers: 73.23.-b, 03.65.Nk, 03.67.Mn 



The scattering matrix, taking asymptotically free in- 
coming states through an interaction region and provid- 
ing the free outgoing states, is of huge basic and prac- 
tical interest. Originally introduced by Born pj] and by 
Wheeler and independently by Heisenberg 0] in atomic 
and particle physics, its application to electron transport 
3 1 has made it into a central tool of mesoscopic physics 
J|. Its formulation for non-interacting electrons pro- 
vides the two-terminal conductance between reservoirs 
in terms of the transmission probability across the scat- 
terer Q and generalizations to complex setups Q have 
been successful. Even more, the scattering matrix ap- 
proach has been used to study noise Q and a description 
of fluctuations in non-equilibrium situations is provided 
by the full counting statistics 01 j with numerous applica- 
tions known today ^ . These successes ask for a general- 
ization to the interacting case for which Goorden et 
al. [lo| have recently derived the scattering matrix of two 
coupled conductors within a perturbative calculation. 

Of generic interest in mesoscopic physics is the quan- 
tum dot, where the transport across is usually studied 
for a finite bias, and different approximate schemes have 
been used to include interactions in the transport anal- 
ysis Uj. Recent interest concentrates on a new type of 
transport, where designed voltage pulses are generated to 
send a finite number of electrons towards the scattering 
region see Ref. for recent experiments in this 
direction. In this situation, novel effects, such as those 
due to particle exchange, show up quite prominently [14 1. 
Moreover, this type of control on individual particles al- 
lows to study fundamental quantum properties such as 
the entanglement between electrons [15| . 

In this letter, we derive the TV-particle scattering ma- 
trix for electrons that propagate freely in the leads but 
exhibit Coulomb repulsion when interacting on the dot. 
Within our formalism, the interaction is accounted for 
by the Hamiltonian Hint = e^7V^/2C, where N is the 
dot's electron number operator and C denotes its capac- 
itance. We use our results to study the wave function 
and degree of entanglement of two scattered electrons. 
Below, we construct the two-particle scattering matrix 
for the quantum dot including Coulomb interaction and 
generalize the result to the iV-particle situation. 



Usually, the scattering matrix connects states at given 
energies; here, we start with the propagator describing 
the scattering of wave packets in coordinate space. We 
start with a (properly symmetrized, spin indices are sup- 
pressed) incident two-electron wave function at time ti, 
^'"(y, tl), with y = {yi,y2}- The scattered wave func- 
tion at later times t2 > ti can be obtained with the help 
of the two-particle propagator K^^\x, 12] y, ti) describing 
the evolution of two particles from the initial positions y 
at time ti to the final positions x at ^2, 

*™*(f,<2)= [ d'yK^'\S,t2;y,h)¥^y,h). (1) 



The two-particle propagator K'-^^ can be defined through 
a Feynman path integral over trajectories x{t), 

K^'\x,t2;y,h)^Jv[x] exp(J ^ 'dt L<^'(f; f )) , (2) 

with the boundary conditions x(ti) = y. Here, U^^{x,x) 
is the system's Lagrangian including kinetic (cx m), dot 
potential (C/), and interaction (cx Uc = 2e^/C) energies, 



[mxf 



U. 



U{xi) --^[xd(a;i) + Xd(a;2)] , (3) 



with the characteristic function Xd{x) of the dot equal to 
unity within the dot and zero outside [l6| . 

Without interaction, the two-particle propagator fac- 
torizes, K'-"'' {x,t2;y,ti) = Il^K^^^{xi,t2;yi,ti) with 
K^^\x,t2]y,ti) the one-particle propagator, while the 
interaction mixes the particle trajectories. A Hubbard- 
Stratonovich transformation [l7| with the real auxiliary 
field z{t) allows us to decouple the quadratic interaction 

K^'\x,t2;y,h) ^ J V[z]exp[i^ J dtz\t)\ (4) 

X ifj'^'l {xi,t2;yi,ti) ifj'^'l {x2,t2;y2,ti), 

where K^^^^{x,t2',y,ti) is the one-particle propagator in 
the presence of a fluctuating potential Uc{t) = Ucz{t), 
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Next, we introduce the scattering matrix S'^^^(e) of the 
dot in the absence of the fluctuating potential Uc{t); the 
indices a,P € {l,r} specify the lead indices for the out- 
going (a) and incoming scattering channels and e de- 
notes the energy variable. We describe the dot through 
the resonance positions (e^) and (identical) widths (F); 
the scattering matrix S'^p{e) then takes the form 

where the constant 2x2 matrices r^^ and s^^ can be 
found from the unitarity conditions. The Fourier trans- 
form provides the real time (t) representation 

3 

where r] = T /Wis the inverse dwell time, utj = ej/Ti is the 
resonance frequency, and 6{t), 9{t) are the usual S- and 
Heaviside functions. The first term in Eq. ^ describes 
the reflection of a particle that has not penetrated into 
the dot, while the subsequent terms correspond to pro- 
cesses where the particle has spent a time r inside the dot; 
the factor e~'"^^'^ describes the accumulated phase. The 
presence of the fluctuating potential Uc{t) contributes an 
additional phase to the one-particle scattering matrix ©, 

= -ti)exp(-^^*' US)dt), (7) 

where ti and ^2 denote the arrival and escape times 
of the particle (we assume escape amplitudes that 



depend weakly on energy). The additional phase 
derives from the gauge transformation (x,t) — > 
^'(a;, i) exp(— I; / dt' Uc{t')); we neglect the scattering 
potential of size ^U^huj 1^ ^ U c\/ vtj ac arising at the 
dot's edge (see Ref. JS; here, and e denote the dot 
size and its dielectric constant, a is the fine structure 
constant). With a typical mesoscopic setup in mind, we 
assume velocities v of order of the Fermi velocity and en- 
ergies larger than the Coulomb energy, e(fc) ^ Uc\ with 
typical ratios v/c ^ 10~^, we can safely ignore the fluc- 
tuation corrections in the dot potential U{x). 



Next, we express the propagator K\])^ through the scat- 
tering matrix ([7]). To simplify matters, we linearize the 
spectrum, e(fc) = fivk] a particle escaped out of the 
dot then never returns. In terms of trajectories, the 
scattering process involves three stages: i) the ballistic 
motion with velocity v towards the dot, ii) the dwell 
time in the dot, and, in) the ballistic propagation away 
from the dot. We define the coordinates in the left 
[x < 0) and right {x > 0) leads with respect to the left 
{x = 0~) and right {x — 0^) dot boundaries and express 
the propagator iiTj^j through the scattering matrix ([7]), 

■^a/3,[z](2;'*2;y,ii) = S''^lxA^T,s)lv, where s = ii + \y\lv 
and T = t2 — \x\lv are the arrival and escape times 
of the particle to and from the dot; similar definitions 
(si = ii -f \yi\lv and Ti =ti — \xi\lv) apply to the two- 
particle scattering matrix, for which we write 



'^aia2/3l/32(^l'^2;Sl,S2) = 5'i\'^^(Ti-Si)S'^'^>^^(T2-S2)(exp|-iWc / dt z{t) + f dt z{t) ^'^ 



(8) 



and in terms of which the two-particle propagator (j4]) assumes the form X*^' — S''^^ /v'^. Here, ojc — Uc/Ti and the 
average in Eq. ^ is taken with respect to the fluctuating Gaussian fleld z{t). The latter is S correlated in time, 
{z{t2)z{ti)) — {i/2LOc)5{t2 — ti) (the complex propagator in Eq. (|4]) generates an imaginary correlator for the real field 
z{t)) and thus the last factor in Eq. ^ can be explicitly averaged over with the result 



•^aias/Sift (n, T2; Sl, S2) = S^l^^ (ti - Si ) 5^ (t2 - S2) exp(-iWcTl2 /2) , 



(9) 



where S\^p is the scattering matrix ((6]) with rcnormalized 
resonance frequencies Cj = ej + C/c/4 and T12 is the time 
the two particles spend together in the dot, 

TI2 = idn - S2I + |t2 - Sl| - |ti - T2\ - \SI - S2I). 

This two-particle scattering matrix ([9]) is the key result 
of this Letter. All effects of Coulomb interaction are ac- 



counted for by rcnormalized resonance energies due to 
self-interaction of individual electrons in the dot and an 
additional phase accumulated by the electrons during 
their simultaneous presence in the quantum dot. 

An inverse Fourier transformation provides us with the 
energy representation 
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The first term describes the non-interacting process 
where particles scatter sequentiahy. The second term 
accounts for inelastic processes where only the total en- 
ergy = £i -I- £2 is conserved. The Coulomb interaction 
generates additional poles at Ejk = e.j + ik + Uc/'2, — iV 
involving the total energy E. These interaction-induced 
singularities cannot be obtained via a perturbative ex- 
pansion for large [/c 3> F. For weak interaction Uc 
or far away from the resonances \E — tj — Uc, the 

expansion of Eq. ([TO|l to first order in Uc reproduces the 
perturbative result obtained in Ref. 

The above derivation for the two-particle scattering 
matrix can be generalized to N particles; the averaging 
over the field z(t) generates an additional phase factor ac- 
counting for the pairwise interaction of particles residing 
simultaneously (for a time Tjk) on the dot, 

N N 

]>k j=i 

The above result also holds true for a multichannel setup, 
with aj, Pj, j = 1,...,N turning into multichannel 
indices. In particular, the results can be straightfor- 
wardly applied to the experimental setup (20l | with two 



parallel leads feeding/emptying two capacitively coupled 
dots that has been recently used to measure interaction- 
induced cross correlations, see also Ref. 

In applying our results to realistic mesoscopic prob- 
lems, we have to avoid mixing between the scattered par- 
ticles and the electrons in the Fermi sea. Hence, we do 
not consider situations with levels within the distance F 
around the Fermi energy e-p and assume that Uc does not 
shift a level across e^; the latter allows us to ignore com- 
plications due to the Kondo effect . In the following, 
we study the scattering problem of two single-electron 
excitations created above the Fermi sea and a quantum 
dot with only one resonance at eg above the Fermi energy 
£f, eo — £f ^ F. The scattering matrix pO)) then tells, 
that (the non-trivial component of) the scattered wave 
function involves energies near eg and e+ = cq + Uc/2. 

We start from a two-electron state with wave function 
^™{xi , X2) created at time i = in the left lead and mov- 
ing towards to the dot [21J . The scattered wave is given 
by Eq. ([1]) and can be expressed in terms of retarded vari- 
ables ^1.2 = |a;i,2| — Vpt, with the Fermi velocity. The 
scattered wave to the right of the dot involving tunneling 
of both electrons assumes the form (Y = yl + j/2) 



+e 



dyi / dy2e~''""e 



dyidy2'^™{vi,y2)e 



ikoY(ix+i2-Y)/i 



-»fce(y-|yi-y2|)/2g-»fcoyg(6-Hf2-i')/f 



;^?(6-6)^'"(yl,2/2) + e(a-6)*'"(y2,2;l) 



(11) 



where £ — 21% /Tv^ is the real-space width of the scattered 
wave, ^> = max{^i,^2}, £,< = min{^i,^2}, fco = wq/wf, 
kc = Wc/2vf, and fc+ = ~\- kc. The second term de- 
scribes the process where the electrons do not overlap in 
the dot, while the term cx e^^+^>e^^°^^ deals with the 
case where both electrons occupy the dot simultaneously 
during scattering. For electrons in a spin-triplet state 
with anti-symmetric orbital wave function ^'"(yi, 1/2), 
this term vanishes and no interaction effects survive, a 
consequence of the Pauli principle. 

Next, we choose a spin-singlet incoming state with one 
orbit 0, hence 4''"(yi, 1/2) = 4'{yi)4'{y2)- To simplify mat- 



ters, we choose an exponentially truncated plane wave 
function with a wave vector in resonance with the dot, 
(t){x) = e{a - a;)e*'=°(^-'')e(^-'')/2L/^^ ^^^^^ a < is 

the initial position of the wave packet at i = and L de- 
notes its width; this choice allows to perform all integrals 
in Eq. pT|) explicitly. We determine the probabilities Pi 
and P2 to transmit one and two electrons through the dot 
(at resonance, where IsrlI = 1) and find them to depend 
only on the two dimensionless parameters a = kcL and 
(3 — £/2L quantifying the interaction and the resonance 
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width, respectively, 



Pi 



P2 



2(3 



1 + 3/3 



1 



2 + 3/3 1 

(l + /3)2 " (a/3)2 + (1 + /3)2 

l + /3(3 + /3) 



(l + 3/3)(l + /3)2 

4(a/3)2+3(/3+3)(/3+l)2 
((a/?)2 + (l+/3)2)(4(a/3)2 + (3 + /3)2) 



,(12) 
(13) 



Without Coulomb interaction (a — 0) or for vanishing 
dwell time (/? = 0), we recover the results Pi — 2p{l ~ 
p) and P2 = p^, with p = 1/(1 + /3) the single-particle 
tunnelling probability. Both, Coulomb interaction (a > 
0) and finite dwell time (/3 > 0) suppress the probability 
P2 as compared with the non- interacting value p^. Even 
infinite Coulomb energy (a — + 00) still permits tunnelling 
of two electrons through the dot via sequential tunneling. 

Finally, we show that the Coulomb interaction in the 
dot leads to an orbital entanglement of the two particles 
(for interaction-induced spin entanglement in a quantum 
dot, see Ref. l22l) . Here, we concentrate on the component 
of the wave function where two electrons are transmit- 
ted to the right and estimate its degree of entanglement, 
which is entirely due to the interaction in the dot. We 
analyze the situation where the length L of the incoming 
wave packet tends to zero, hence (3 ^ 00. In this case the 
normalized wave function on the right has the universal 
form (independent of ^f'") 



where ^1^2 = 2:1,2 — a— fpi < 0. Eq. describes a two- 
electron state with different momenta fc+ and fep < fc+, 
as has to be expected since the first electron escaping 
carries an energy shifted up by the Coulomb interaction. 
The state (1141) can be rewritten in a form 



(15) 



reminding about the original Einstein-Podolsky-Rosen 
state with orbital entanglement [1^. To quantify its 
entanglement, one may calculate the von Neumann 
entropy E of the reduced density matrix p{x,x') = 
J da;2 ^'rr(x, X2)5'Rpt(a;', a;2). Instead, we determine the 
purity Tl{p) = trp^, which is unity for separable states 
and provides the lower limit E > 1 — 11. With A = 
ikc£/{2 — ikci), we find the density matrix 

p{^,a = (2/^)(?(-0^(-e')e«+«''/'e^'="(«-«') (16) 
X [1 + A9{^-aie^i^' - e^'==(«-«')e2«'/«) 
+A*(?(e'-0(e'«'/' - e*'==«-«')e2«/^)] , 

that results into a purity H = [1 + 2/(1 + {kJ/4)'^)]/3. 
We conclude that at finite Uc the state (fTSj) is entan- 
gled and the degree of entanglement saturates as the 
Coulomb interaction becomes larger than the resonance 



width, kci = Uc/T ^ 1, i.e., when the energies of the 
escaped particles become distinguishable. 

In conclusion, our expression for the multi-particle 
scattering matrix accounts for the redistribution of par- 
ticle energies during inelastic scattering with the appear- 
ance of new resonance poles that cannot be obtained per- 
turbatively for large Uc ^ F. As an application, we 
have studied the case where two electrons are transmit- 
ted across a dot with a single resonance and have investi- 
gated the ensuing orbital entanglement and the reduction 
in the two-particle transmission due to the interaction. 

We thank Fabian Hassler and Konstantin Matveev 
for discussions and acknowledge the financial support 
from the Swiss National Foundation (through the pro- 
gram MaNEP and the CTS-ETHZ), the Russian Founda- 
tion for Basic Research (08-02-00767a), and the program 
'Quantum Macrophysics' of the RAS. 



[1] 
[2] 

[3] 

[4] 



[5] 

*rr(6,6) = (2/^)e''=+«>e''^'°«< e(«l+«^)/^ (14) [6] 

[7] 



[8] 
[9] 

[10] 
[11] 



[12] 



[13] 
[14] 

[15] 



[16] 



M. Born, Z. Physik 37, 863 (1926). 

J.A. Wheeler, Phys. Rev. 52, 1107 (1937); W. Heisen- 

berg, Z. Physik 120, 513 and 673 (1943). 

R. Landauer, IBM J. Res. Dev. 1, 223 (1957); R. Lan- 

dauer, Philosoph. Mag. 21, 863 (1970). 

Y. Imry and R. Landauer, Rev. Mod. Phys. 71, S306 

(1999) . 

M. Bxittiker, Y. Imry, R. Landauer, and S. Pinhas, Phys. 
Rev. B 31, 6207 (1985). 

G.B. Lesovik, JETP Lett. 49, 592 (1989); M. Biittiker, 
Phys. Rev. Lett. 65, 2901 (1990). 

L.S. Levitov and G.B. Lesovik, JETP Lett. 58, 230 
(1993); L.S. Levitov, H. Lee, and G.B. Lesovik, J. Math. 
Phys. 37, 4845 (1996). 

Ya.M. Blanter and M. Biittiker, Phys. Rep. 336, 1 

(2000) ; Noise m Mesoscopic Physics, Vol. 97 of NATO 
Science Series, ed. Y. Nazarov (Kluwer, Delft, 2003). 

A. Aharony, O. Entin-Wohlman, Y. Levinson, and Y. 
Imry, Ann. Phys. (Leipzig) 8, 685 (1999); A. Aharony, 
O. Entin-Wohlman, and Y. Imry, Phys. Rev. B 61, 5452 
(2000). 

M.C. Goorden and M. Biittiker, Phys. Rev. Lett. 99, 
146801 (2007). 

M. Kindermann and Yu.V. Nazarov, Phys. Rev. Lett. 
91, 136802 (2003); M. Kindermann and B. Trauzettel, 
Phys. Rev. Lett. 94, 166803 (2005); A. Komnik and A.O. 
Gogolin, Phys. Rev. Lett. 94, 216601 (2005); Y. Utsumi, 
D.S. Golubev, and G. Schon, Phys. Rev. Lett. 96, 086803 
(2006). 

D.A. Ivanov and L.S. Levitov, JETP Lett. 58, 461 (1993); 
J. Keeling, I. Klich, and L.S. Levitov, Phys. Rev. Lett. 
97, 116403 (2006). 

G. Feve et al, Science 316, 1169 (2007). 

F. Hassler, G.B. Lesovik, and G. Blatter, Phys. Rev. Lett. 

99, 076804 (2007). 

A.V. Lebedev, G.B. Lesovik, and G. Blatter, Phys. Rev. 
B 72, 245314 (2005); P. Samuelsson and M. Biittiker, 
Phys. Rev. B 71, 245317 (2005); C.W.J. Beenakker, M. 
Titov, and B. Trauzettel, New J. Phys. 7, 186 (2005). 
The capacitive form implies an instantaneous and uni- 



5 



form change of the dot potential upon tunnehng of the 
electron; the approximation is valid for dots with plasma 
frequencies well above ballistic ones (of order GHz). 

[17] J.W. Negele and H. Orland, Quantum Many-Particle 
Systems, Frontiers in Physics Lecture Note Series (Ad- 
dison Wesley, Redwood City CA, 1988), p. 333. 

[18] Eq. ([4]) gives typical fluctuations ~ h/Ucdt, with bal- 
listic propagation giving the fastest time 5t ~ Id/y- 

[19] L.I. Glazman and M.E. Raikh, JETP Lett. 47, 452 



(1988). 

[20] D.T. McClure et al, Phys. Rev. Lett. 98, 056801 (2007). 
[21] Use of unit-flux voltage pulses guarantees that excited 

particles are not entangled with the Fermi sea. 
[22] W.D. Oliver, F. Yamaguchi, and Y. Yamamoto, Phys. 

Rev. Lett. 88, 037901 (2002). 
[23] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 

777 (1935). 



